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I.  Introduction 


— ^^The  objective  of  the  present  investigation  is  to  make  a  comprehen¬ 
sive  experimental  evaluation  of  a  recent  theory  of  elastic-plastic 
impact^^  and  to  modify  the  theory  to  provide  a  more  realistic  description 

of  the  impact  characteristics.  Previous  work  has  included  investigation 

V  f  j  - 

of  subsonic  localized  impact  damage  in  glass  ,  ZnS  ’  ’ '  ,  transforma- 

V-l  f  1  .11}— ~ — i  /r> - 

tion  toughened  zirconia  (TT-2r02? ’  ~ T  ,  and  other  ceramics 


A  major  advantage  of  the  theory  is  that  it  provides  a  means  to  predict  the 
response  of  a  material  subjected  to  impact  loading  using  contact  character¬ 
istics  determined  by  static  loading  tests.  The  assumption  that  the  damage 
induced  at  a  particular  maximum  static  load  is  the  same  as  that  resulting 
from  an  equal  maximum  impact  load  has  frequently  been  invoked  or  su/^gr  ct- 
ed^^  ’■*■**)  ^  A  recent  report  describes  similarities  and  differences  in 
damage  induced  by  static  and  impact  loading^^T"  A  preliminary  evaluation 
of  the  theory  of  elastic-plastic  impact  using  limited  data  showed  that  the 
theory  yields  promising  results  ' .  The  present  report  describes  certain 
modifications  of  the  theory  and  a  more  complete  evaluation  of  the  theory 
based  on  more  extensive  data.* 


II.  Experimental  Procedures 

Glass,  CVD  ZnS ,  and  TT-ZrO^  plates  described  in  Table  I  were  pre¬ 
pared  as  follows: 

Glass — Larger  glass  plates  were  cut  using  a  glass  cutter  to  form 
the  square  specimens. 

ZnS — The  square  plates  were  used  in  the  polished  condition  in  which 
they  were  secured  from  the  manufacturer. 

TT-ZrO^ — The  calcia  partially  stabilized  zirconia  plates  were  pre¬ 
pared  at  CSIRO  where  they  were  aged  at  1300°C  for  various  periods  of  time 
to  form  three  grades  in  which  the  K  increased  with  increasing  aging  time 

1  v 

(as  fired,  0  hours;  under  aged,  30  hours;  peak  aged,  48  hours).  At 
Ceramic  Finishing  Company  the  TT-ZrC^  plates  were  polished  on  one  surface 
using  1/4  pm  diamond  powder  as  the  last  step.  These  specimens  were  frac¬ 
tured  in  earlier  experiments  and  the  remaining  ends  of  these  specimens  were 
used  in  the  present  experiments. 

The  physical  properties  of  the  materials  are  listed  in  Table  II. 

The  plates  were  coated  with  a  thin  layer  of  soot  by  moving  the  plates 
back  and  forth  through  a  candle  flame  until  a  layer  of  sufficient  thickness 
was  formed.  The  imprint  of  the  sphere  in  the  soot  layer  was  used  to  measure 
the  contact  radius.  The  assumption  that  the  radius  of  the  imprint  repre¬ 
sents  the  contact  radius,  including  the  effects  of  variation  in  soot  layer 
thickness,  is  discussed  in  a  later  section. 


Table  I.  Description  of  Specimens. 


TabLe  II.  Material  Properties. 


calculations  involving  equations  for  the  elastic  case  (see  Table  III), 
calculations  involving  equations  for  the  elastic  case  (see  Table  III). 
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The  specimens  were  subjected  to  either  static  or  impact  loading.  The 
static  loads  were  applied  using  a  table  top,  Instron  testing  machine. 

When  the  glass  spheres  were  used  as  indenters,  the  spheres  were  fastened 
to  the  crosshead  and  the  glass  or  ZnS  plates  were  placed  on  the  load  cell. 
The  tungsten  carbide  spheres  were  fastened  to  the  load  cell  and  the  TT-ZrO^ 
specimens  were  fastened  to  the  cross  head  because  of  problems  involved  in 
gluing  the  very  small  spheres  to  the  cross  head. 

For  impact  loading,  the  spheres  were  coated  with  a  thin  layer  of  alum¬ 
inum  to  improve  their  reflectivity  and  accelerated  using  a  gas  gun  with 
compressed  gas.  The  specimens  were  mounted  on  a  steel  plate  with  a 
large  surface  perpendicular  to  the  oath  of  the  sphere.  The  impacting  and 
rebounding  velocities  of  the  spheres  were  measured  using  photographs  taken 
using  a  camera  with  anopen  shutter.  The  spheres  were  illuminated  at  known 
time  intervals  using  a  stroboscope  and  the  velocities  were  calculated  by 
dividing  the  distance  between  images  by  the  time  interval . 

The  contact  and  indentation radi i  were  measured  using  a  grid  in  the 
eyepiece  of  a  microscope. 
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III.  Review  of  the  Theory 

n 

The  theory  of  elastic-plastic  impact is  based  on  the  following 
empirical  relations  for  static  loading: 

n  =  *(r  /r  )C  (1) 

c  o 

n'  =  k*  (r./r  )'•'  (2) 

1  o 

where  n  and  n'  are  the  contact  and  indentation  pressures,  respectively; 

r  is  the  contact  radius,  r  is  the  indentation  radius,  and  r  ,  >  '  ,  •'  and 
c  1 

are  contact  characteristics  that  are  determined  from  log- log  plots  of 

n  vs  r  /r  and  n'  vs  r./r  . 
co  l  o 

The  characteristics  f,  and  f, '  account  for  work  hardening  and  other 

material  characteristics  that  control  the  deformation  under  static  loading 

conditions.  The  characteristic  k  represents  the  pressure  at  which  the 

contact  radius  equals  the  radius  of  the  spherical  indenter.  Likewise,  >; ' 

represents  the  pressure  at  which  the  indentation  radius  equals  the  radius 

of  the  spherical  intenter.  At  high  r/r  ratios,  Equations  (1)  and  (2)  no 

o 

longer  apply  because  of  deformation  of  the  sphere  and  increasing  errors 
in  the  geometrical  assumptions.  Also,  for  a  given  material  we  would 
expect  k  to  he  less  titan  because  r.  is  less  than  r  for  all  r  <  r 

tco 

so  that  the  indentation  area  is  less  than  the  contact  area  and  n  is  less 


than  n ' . 
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The  pressures  are  taken  to  be  the  load  divided  by  the  projection  of 
the  area,  so  that  (1)  and  (2)  can  also  be  written  in  the  form 

i  r  + ) 

P  =  ner  "(r  /r  )  '  (3) 

O  CO 

2  r  •  +9 

l1  =  nK'r  (r./r  )tj  (4) 

o  l  o 

An  expression  for  the  impact  loading  contact  radius  can  be  derived  by 
assuming  that  the  sphere  remains  rigid.  Then,  the  work  done  by  the  load 
as  given  in  Equation  (3),  to  the  point  of  maximum  penetration,  can  be 
equated  to  the  initial  kinetic  energy  of  the  impacting  sphere.  This 
y ields 


r 

c 


r 

o 


M(E+4) 

.  3 

dsn  r 

o 


X  2 

E+4  r,+4- 

v 

o 


(5) 


By  assuming  that  equal  static  and  impact  loads  result  in  equal  contact 

radii  we  can  substitute  for  r  in  liquation  (3)  using  Equation  (3)  to  get 

c 

the  following  expression  for  the  impact  load 


I’ 


2 .  M  (  f,+4  ) 

ro  [7 - 3 

2k  it  r 


f.+2  2  (f.+2  ) 

E+4  (f.+4) 

V 

o 


(6) 
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Similar  manipulations  provide  expressions  for  t lie  indentation  radius, 

penetration  time,  depth  oi  damage,  remaining  strength,  contact  time,  and 

coefficient  of  restitution.  The  equations  used  in  the  present  work  are 

given  in  Table  111  which  also  includes  equations  for  completely  elastic 

(17  ?7) 

response  which  will  be  referred  to  as  tlu-  Hertzian  case  ’ 

Modification  of  the  Analyses 

Roesler  and  others^^’^^  used  dimensional  analysis  and  energy  argu¬ 
ments  to  derive  equations  describing  tile  variation  of  crack  size  (c)  due 
to  varying  contact  load  (P)  in  the  "far  field"  where  the  effects  of  the 

nature  of  the  particular  contact  arc  considered  to  be  small.  In  the  form 

(20) 

used  by  Lawn  and  Marshall  this  equation  in 


P 


Kic  3/2 
- c 


(7) 


in  which  Kr  is  the  critical  stress  intensify  factor  and  v  is  a  constant 
Lc 

that  accounts  for  geometrical  and  :Y:et  io::-i!  ep.-n . .  riot  i.\-  of  :  Ijt  ■  : 

The  equation  was  incorporated  into  the  original  analysis  of  elastic-plastic 
impact  to  permit  calculation  of  the  depth  of  damage  and  the  remaining 

strength.  However,  in  an  investigation  of  penetration  of  damage  during 

,  ,  .  (21) 

single  point  diamond  grinding  ot  ceramics,  K i rchne r ,  (iruvor,  and  Richard 

found  that  for  diamond  points  with  varying  sharpness,  the  depth  of  d. image 

(c)  varied  systematically  with  the  length  of  contact  even  in  the  "far  field' 

where  the  effect  of  contact  characteristics  was  expected  to  he  minimal. 


List  of  Svmbols 


Contact  pressure 

Indentation  pressure 

Resistance  to  Penetration 

Resistance  to  indentation 

Modified  < 

Exponent  characteristic  of  the  variation  of  the  contact 
pressure  with  increasing  load 

Exponent  related  to  work  hardneing 

Modified  f, 

Radius  of  impacting  sphere 
Mass  of  impacting  sphere 

Contact  radius 

Indentation  radius 

Maximum  load  perpendicular  to  the  surface 

..  ,  ,  ,  i  -  1  for  sphere 

Voung  s  moduLus  .  _  _ 

]  =  2  for  plate 

Poisson's  ratio  j  =  i,  2  as  above 
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K  . 


a,  6  - 


Ic 


Q.-i-Kl-vp  + 


Impact  velocity 


(1-  v.) 

K  .  =  —  re  — 

1 


Coefficient  of  restitution 
£+2 
£+4 


ft  -  <ur 


2.M(C+4) 


°  2<nr  3 
o 


Indentation  diameter 


Crack  length 


Constants  determined  from  static  data 
where  P/d.  =  a  + 


Critical  stress  intensity  factor 


Proportionality  factor  accounting  for  unknown  geometrical  and 
frictional  factors 


Depth  of  the  disturbed  zone 
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Based  on  L ho  above  observations,  Conway  and  Kirchner  analyzed 

penetration  of  surface  damage  for  the  case  in  which  a  "disturbed"  zone 
is  formed  under  the  contact.  This  "disturbed"  zone  may  be  formed  as  a 
result  of  plastic  deformation,  shear  cracking  or  crushing  under  the  con¬ 
tact.  This  analysis  yielded  a  relation  of  the  form 


P  “  z  c1^2  (8) 


in  which  z^  is  the  depth  of  the  disturbed  zone.  A  proportionality  factor 

can  be  introduced  which  accounts  for  other  characteristics  of  the  contact 

and  variations  in  KT  from  one  material  to  another.  Similar  relations 
Tc 

might  arise  as  a  result  of  line  contact  loading,  wedging,  or  residual 
stresses . 

Shear  failures  under  contacts  have  characteristic  flow  or  crack 

(23) 

patterns  that  follow  shear  stress  trajectories  .  Based  on  the  shapes 

of  these  patterns  there  appears  to  be  a  direct  relationship  between  the 

depth  of  the  disturbed  z.one  and  the  contact  width.  Therefore,  it  may  be 

appropriate  in  many  cases,  such  as  single  point  diamond  grinding  with 

diamonds  having  various  degrees  of  sharpness ,  to  substitute  the  contact 

widtli  for  z  in  the  above  relation, 
o 

A  "disturbed"  zone  has  been  observed  in  ZnS  statically  indented  by 

glass  spheres.  When  the  indentation  diameter  (d  )  was  substituted  for  z^ 

1/2 

in  the  above  relation  and  I’/d^  was  plotted  vs  c  a  good  linear  fit  was 
(8,9) 


observed 


Based  on  the  above  experience,  relation  (8)  has  been  suits t i tuted  for 
Equation  (7)  in  the  original  derivations  of  the  equations  for  calculating 
the  depth  of  damage  and  the  remaining  strength  after  impact. 

Remedies  for  a  Problem  in  Measuring  r 
_ c^ 

As  mentioned  previously,  one  should  expect  <  to  be  less  than  <’  because 
r^_  >  r^  in  all  cases.  However,  experimentally  one  finds  that,  as  shown  in 
Figure  1,  the  extension  of  the  data  representing  the  contact  tends  to 
intersect  the  data  representing  the  indentation  at  a  point  r  <  r^  rather 
than  at  r  =  r^.  Two  possible  explanations  of  the  observation  were  con¬ 
sidered: 

1.  The  elastic  recovery  in  the  plane  of  the  specimen  surface  may 

( ?8) 

decrease  r^  substantially 

2.  The  increased  discontinuity  at  the  rim  of  the  indentation  with 

increasing  indentation  size  may  cause  r^  to  increase  less  with 

load  than  otherwise  expected.  In  terms  of  the  graph  this  implies 

a  shallower  slope  than  the  original  contact  radius  curve  and  an 

intersection  at  r  5  r  . 

c 

It  is  unfortunate  that  the  cause  of  the  difficulty  is  not  understood  because 
the  remedies  implied  bv  the  above  explanations  are  quite  different.  In  the 
first  case  it  is  reasonable  to  assume,  because  the  yield  stress  is  ex¬ 
ceeded,  that  the  correct  curve  is  simply  shifted  to  the  right  as  indicated 
in  Figure  2,  so  that  it  intersects  at  r  =  r  .  This  is  the  same  as  substi¬ 
tuting  k'  for  K .  In  the  second  case  the  data  at  low  r  /r  would  be  considered 

c  o 

to  be  relatively  good  but  the  errors  would  become  greater  at  higher  values 


Pressure  -  GPa  c  Pressure -GPa 


of  r  /r  .  Therefore,  in  this  ease  it  is  reasonable  to  connect  Liu*  good  data 
c  o 

points  at  Low  r  /r  values  witli  r  =  r  as  shown  in  Figure  1. 
o  o  e  o 

In  any  event  it  was  decided  to  assume  that  ^  =  x'  and  to  refer  to  this 
case  as  the  modified  elastic-plastic  theory.  One  would  expect  that  if  the 
other  choice  had  been  made  the  results  would  have  fallen,  in  most  cases, 
between  those  of  the  original  and  modified  theories. 


Pressure  -  GPa 


•l'  V  •  .  r'lfV-  ,  'I  V:  .  !" .  /  r  ,  •('.  I  !  '  i  I  '| 

r  Jr  l  vhi-h  'i:  ID'  Uki'  if  i :  ■  i  1 1  *'  •  !  i  :••••  iii*  i  f :  J  i  1  i  • 

i  ft' 1  ■!  •  i  IT  '  •  I '  I  >  ’ll  I  'l !  i  '  >ll  fill  i  i  r'lll::.  :i  1  ■  ■  r  t  i.-ili  XI’ 

iMti'.’l'  •  V  f  it)  i  ’If.  . 


IV.  Determination  of  c,  k'  ,  C ,  and  f, ' 


1  > 


In  the  case  of  the  r-tut.  i  •••• !  loaded  r.o- >:•  i ;•>  :  ,  n  'nn  n'  d«  *  •  •  rir.  i  ? .• 

by  dividing  the  applied  load  by  the  contact  area  and  the  indentation  area, 

respectively.  Then,  n  was  plotted  vs  (r  /r  )  and  n'  was  plotted  vs  r./r 

co  to 

on  log-log  plots  as  suggested  by  equations  (1)  and  (2).  K  ,  f, ,  and  i,' 

were  calculated  from  the  slopes  of  these  plots  (Table  IV). 

The  original  data  were  less  accurate  than  desired.  The  principal 

problem  involved  measurement  of  the  contact  radius.  The  soot  layer  which 

is  used  to  measure  r^  must  be  thick  enough  to  provide  a  distinct  contrast 

at  the  edge  of  the  impression  so  that  an  accurate  measurement  can  be  made, 

but,  if  the  soot  layer  is  too  thick,  the  sphere  will  pick  up  soot  from 

outside  the  contact  area.  In  some  cases  this  results  in  an  irregular 

boundary.  The  minimum  radii  were  measured  in  these  cases.  Because  the 

soot  layer  was  applied  manually  the  thickness  varied  to  some  extent  from 

specimen  to  specimen.  The  data  points  representing  eases  in  which  t  tic 

soot  layer  is  definitely  too  thick  are  noted  on  some  ot  the  graphs.  In 

general,  the  r,  measurements  are  more  accurate  and  less  scattered  than 

the  r  measurements.  One  benefit  of  substituting  >  '  tor  r  as  is  done  in 
c 

the  modified  theory  is  that  it  suppresses  some  of  the  effect  ot  inaccuracies 
in  the  r  measurements. 


2  i 


V.  Experimental  Evaluation  of  the  Elastic-Plastic  Impact  Theory 

Class  Spheres  on  Glass  Plates 

This  material  combination  involving  glass  spheres  and  glass  plates 
represents  an  end  member  case  in  which  the  response  at  low  impact  velocities 
is  almost  completely  elastic.  Therefore,  one  would  expect  the  Hertzian 
elasticity  theory  to  provide  a  good  description  of  the  impact.  The  contact 
characteristics  were  £  =  2.2  and  k  =  126  GPa.  The  glass  does  not  indent 
permanently  so  £ '  and  <'  were  not  determined. 

The  contact  radii  predicted  for  various  impact  velocities  are  compared 
with  the  Hertzian  theory  predictions  in  Figure  4.  The  Hertzian  theory 
fits  the  data  best  as  expected  but  both  theories  yield  reasonable  approxi¬ 
mations.  At  the  highest  impact  velocities  above  100  ms  '  crushing  occurred, 
increasing  the  contact  radii  above  the  values  otherwise  expected.  At  inter¬ 
mediate  velocities  ring  and  cone  cracks  formed  in  many  specimens.  There  is 
no  definite  evidence  that  the  presence  of  these  cracks  affected  the  contact 
radii . 

The  loads  predicted  for  various  impact  velocities  are  plotted  in  Figure 
6.  The  predictions  of  the  two  theories  agree  at  low  velocities  but  at 
higher  velocities  the  elastic-plastic  prediction  leads  to  higher  loads 
than  the  Hertzian  theory.  Again,  because  the  response  of  the  glass  is 
elastic,  one  would  expect  the  Hertzian  theory  to  yield  the  best  results. 

The  impact  indentation  radii  and  coefficients  of  restitution  were 
not  calculated  for  impacts  of  glass  spheres  on  glass  plates  because  the 


absence  of  indentation  prevented  determination  of  k'  and  f,'. 


+  Earlier  data 
R.  M.  Gruver,  p4l 
April,  1976  Notebook 
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Glass  Spheres  on  /.nS  Plates 


In  contrast  to  the  case  of  glass  on  glass,  impacts  of  the  glass 
spheres  on  the  ZnS  plates  yield  an  elastic-plastic  response.  Under 
static  loading  conditions,  indentations  were  first  observed  at  a  40  N  load 
but  a  load  of  72  N  was  necessary  to  obtain  a  measurable  indentation.  Inden 
tations  were  observed  at  all  impact  sites.  Impact  velocities  less  than 
about  19  ms  ^  were  impractical  because  the  curvature  of  the  trajectory 
made  it  difficult  to  hit  the  target. 

The  static  contact  characteristics  are  given  in  Table  IV.  The  recent  1 
reported  values  which  were  =  0.71,  k=  4.9  GPa,  f, '  =  0.14,  and  <  '  = 

3.0  GPa  were  used  to  calculate  various  predicted  curves  which  are  compared 
with  the  Hertzian  theory  and  experimental  results.  In  some  cases  *  was 
assumed  to  he  equal  to  < '  as  described  in  Section  111.  These  results 
are  called  the  modified  elastic-plastic  theory  results. 

The  contact  radius  results  plotted  in  Figure  6  show  that  the  contact 
radii  are  only  slightly  larger  than  the  Hertz  theorv  predictions  and  sub¬ 
stantially  less  than  the  other  predictions.  There  are  several  factors 
contributing  to  errors  in  these  predictions  including  the  strain  rate 
dependence  of  the  flow  stress,  neglect  of  the  stress  wave  energv,  neglect 
of  the  fracture  energy,  and  so  forth.  The  above  data  supports  our  earlier 
conclusion  that  the  most  important  source  of  error  is  the  strain  rate 
dependence  (if  the  flow  stress'1  .  As  the  strain  rale  increases,  the  1  low 
stress  increases  and  the  material  responds  more  elasticallv  than  it  does 


under  the  static  loading  conditions  used  to  determine  the  contact  character¬ 


istics. 


This  more  elastic  response  shifts  the  results  toward  the  Hertzian 


600 


o  Experimental  data 


f 

Figure  6  Contact  radius  vs.  impact  velocity  (ZnS,  1  mm  diameter  glass 
s  plieros  )  . 
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theory  curve.  Lt  is  interesting  that  this  shift  occurs  so  uniformly 
over  the  range  of  velocities.  In  this  connection,  it  is  well  to  remember 
that  each  impact  results  in  a  wide  range  of  strain  rates.  During  the  early 
stages  of  the  impact  process  the  highest  strain  rates  are  encountered  but 
these  high  strain  rates  involve  a  relatively  small  volume  of  material. 

Later  in  the  process,  the  strain  rates  are  lower  and  involve  larger  volumes 
of  material.  Apparently,  the  effect  of  averaging  over  these  various  strain 
rates  and  volumes  of  material  is  a  smaller  impact  velocity  dependence  of 
the  contact  radius  than  might  otherwise  have  been  expected. 

The  impact  indentation  radii  calculated  using  tin?  elastic-plastic 
theory  and  the  modified  elastic-plastic  theory  are  plotted  in  Figure  7 
where  the  results  are  compared  with  the  experimental  data.  The  calcula¬ 
tions  overestimate  the  indentation  radii  mainly  because  the  strain  rate 
dependence  of  the  flow  stress  increases  the  elastic  response  compared  with 
that  expected  based  on  the  static  contact  characteristics. 

The  predicted  impact  loads  are  plotted  in  Figure  8.  The  impact  load 
was  not  measured  so  there  are  no  experimental  data  for  comparison.  There¬ 
fore,  comparison  is  made  primarily  with  the  Hertzian  theory.  In  the 
elastic-plastic  case,  yielding  occurs  at  the  contact  which  decreases  the 
average  pressure  in  the  contact  compared  with  that  expected  during  a  com¬ 
pletely  elastic  response.  Yielding  also  increases  the  penetration  of  the 
sphere  into  the  target.  The  net  result  is  that  the  maximum  impact  force 
or  load  is  less  in  the  elastic-plastic  case  than  it  would  he  in  the  elastic 
case.  The  values  indicated  for  the  elastic-plastic  theory  seem  to  be  too 
high  in  relation  to  the  Here theory  hut  those  for  the  modified  elastic 
plastic  theory  seem  reasonable  enough  although  we  do  not,  as  yet,  have  an 


independent  means  of  evaluating,  them. 
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One  possible  explanation  of  t he  high  calculated  values  of  impact 
load  is  that  the  theory  overestimates  the  contact  radii  (Figure  6)  so 
that  the  contact  areas  are  overestimated.  Wien  these  areas  are  multiplied 
by  the  contact  pressures,  the  resulting  load  values  are  overestimated. 

In  an  attempt  to  evaluate  this  effect  the  experimental  contact  radii 
were  combined  with  the  contact  pressures  to  calculate  the  impact  loads 
with  the  results  shown  in  Figure  9.  The  resulting  loads  are  about  50% 
as  great  as  those  represented  by  the  theoretical  curves.  While  this  result 
seems  reasonable  there  is,  so  far,  no  independent  way  to  check  this 
resul t . 

For  impacting  spheres,  the  coefficient  of  restitution  measurements 
provide  a  sensitive  means  to  assess  the  energy  losses  that  occur  during 
the  impact  process.  In  ideal  elastic  impact,  the  sphere  rebounds  at  the 
original  impact  velocity,  the  coefficient  of  restitution  is  one,  and  there 
are  no  energy  losses.  In  an  actual  elastic  impact,  there  are.  stress  wave 
energy  losses,  the  sphere  rebounds  at  a  velocity  slightly  less  than  the 
impact  velocity  and  the  coefficient  of  restitution  is  slightly  less  than 
1,  frequently  in  the  range  from  0.9  to  one.  In  these  cases,  the  measure¬ 
ment  of  the  coefficient  of  restitution  is  a  relatively  sensitive  means 
of  estimating  the  stress  wave  energy.  When  there  is  an  elastic-plastic 
response  the  situation  becomes  more  complicated  because  both  indentation 
energy  and  stress  wave  energy  contribute  to  the  losses  and  tiie  coefficient 
of  restitution  may  decline  to  very  low  values.  Nevertheless,  because  of 
the  sensitivity  of  the  velocity  in  responding,  to  the  energy  losses,  coef¬ 
ficient  of  restitution  measurement s  represent  an  effective  means  to  determine 
how  we  1  1  a  model  corresponds  to  real  itv. 


:j  i 


The  coefficients  of  restitution  calculated  using  the  elastic-plastic 
theory  and  the  modified  theory  at  various  velocities  are  plotted  in  Figure 
10  where  they  are  compared  with  the  experimental  data.  The  predicted  curves 
greatly  underestimate  the  measured  values  apparently  because  the  strain  rate 
dependence  of  the  flow  stress  causes  the  material  to  respond  much  more 
elastically  than  would  be  predicted  based  on  the  response  under  static 
condit ions . 

One  of  the  principal  objectives  of  the  present  program  is  to  predict 
the  crack  sizes  and  the  remaining  strength  after  impact  damage.  Figure  11, 
which  is  based  on  the  relation  P/d;  <*  c37^,  was  presented  in  a  recent 

(  g) 

report  .  The  data  for  c  >  d^,  where  d  is  the  contact  diameter  which  is 
used  to  define  the  boundaries  of  the  near  field,  were  fitted  by 

P/d j  =  2.74  x  107  c1/2  +  7.87  x  K)3  (9) 


Substituting  the  above  equation  in 


P/d  j 
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yields  an  liquation  that  can  he  used  to  predict  radial  crack  length.  The 
results  of  these  theoretical  calculations  are  compared  with  the  experimental 
data  in  Figure  12.  Although  the  modified  theory  yields  an  improvement  in 


the  fit  to  the  experimental  data,  it  is  clear  that  the  theory  still  dot's 


°  °o° 


»  0.60 

Q) 

q: 

h-  0.50 

o 

-  0.40 

c 

<u 

z  0.30 
?  0.20 


Modified  elostic  -  plastic  theory 
Elostic  -  plastic  theory 


0  10  20  30  40  50  60  70 

Impact  Velocity-m-s 


rii'iir*'  10  Tli-  i  vv  Tf  i  <•  i  > -nt,  of4  r<  ■ : ;  t .  i  I  ut.  i  on  v.4..  i.l  ■  import 
i  ni  |  i  •  1.  ■  ■  t  i  !■;/  -1  mm  il  i  am*  ‘L*  t4  r  l-4  :4pir 


lor  it,;/.  ( I4,! 


/■  i  fu  r-  ■  I )•»!  <•  •!'  1  <•»>»:!  b  vr.  i  !!ji •  *  •/•■>•  •  i  1  v  r  '  1  •  :  •  - :  '  . 

'.i:'  i  t.!i* -  rn>  >  I  i  !’ i  -  i  •  ■  ;  '■  i..  ■  ,  .  it  ; 

(l|)  -in  i  (10).  (/'ik';  -l.U  itiiii  ■.]  i  .’jv  '•  i*  -f  I  :ii'  i-  !  . 


35 


not  account  for  the  strong  variation  in  maximum  radial  crack  length  with 
impact  velocity  represented  by  these  data. 

Tungsten  Carbide  Spheres  on  TT-ZrO^ 

Impacts  of  tungsten  carbide  spheres  on  TT-ZrO,;  yield  an  elastic- 
plastic  response  despite  the  fact  that  zirconia  is  much  harder  than  zinc 
sulfide.  The  log  contact  pressure  vs  log  normalized  contact  radius  curves 
for  the  three  grades  of  TT-ZrO^  show  leveling  off  of  the  contact  pressure 
at  normalized  contact  radius  values  above  0.21.  This  behavior  was  not 
observed  in  the  other  materials.  Also,  it  was  not  observed  in  the  inden¬ 
tation  radius  measurements  on  TT-ZrO  1? .  The  mechanism  responsible  for  the 
behavior  is  not  understood.  The  possibility  that  it  might  occur  as  a  result 
of  plastic  deformation  of  the  sphere  was  considered.  The  spheres  wore 
examined  by  optical  microscopy  and  no  evidence  of  plastic  deformation  was 
observed.  However,  it  is  quite  possible  that  small  plastic  deformations 
could  have  gone  undetected.  Another  possibility  is  that  it  might  he  caused 
!>v  <i  phase  transformation,  perhaps  inonoclinic  to  tetragonal  under  the  very 
h  i  v;h  pressure  conditions  existing,  under  the  contact.  With  regard  to  both 
of  these  mechanisms,  it  is  not  elear  whv  the  leveling,  off  is  observed  only 
in  the  contact  pressure  and  not  in  the  indentation  pressure. 

The  contaet  i ha rac t er i s t i cs  measured  for  the  three  grades  of  TT-ZrO 
can  he  compared  (Table  IV).  Although  there  seem  to  he  systematic  variations, 
in  view  ol  the  unee rl a i n t ies  in  the  results,  these  variations  may  not  he 
significant.  The  indentation  radii  observed  .it  various  static  loads  on 
as-fired  and  peak-aged  TT-ZrO  are  compared  in  figure  )  1.  I.itt  le  difference 


) 


was  observed  in  the  results  for  the  two  materials. 


Therefore,  the  aging 


Load 
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treatment  does  not  seem  to  have  an  appreciable  effect  on  the  plastic 
flow  properties  of  TT-ZrO^. 

The  contact  radius  vs  impact  velocity  predictions  are  plotted  in 
Figures  14-16  where  they  are  compared  with  experimental  data.  Because 
the  measured  contact  radii  increase  w'i  th  the  thickness  of  the  soot  layer 
but  the  theory  is  based  on  the  actual  contact  witli  the  specimen  surface, 
it  seems  reasonable  to  give  added  weight  to  the  lower  contact  radius  values 
in  interpreting  these  results.  Using  this  t;  prong’n ,  *  !cj  orii'tvil  and 
modified  elastic-plastic  theories  overestimate  the  contact  radii.  Again, 
it  seems  likely  that  the  contact  radii  arc  smaller  than  expected  because 
the  strain  rate  dependence  of  the  flow  stress  causes  the  T'i-ZrO.,  to  respond 
more  elastically  than  expected  based  on  the  contact  parameters  determined 
during  static  loading.  It  is  also  possible  that  the  data  are  slightly 
underestimated  as  a  result  of  the  elastic  contraction  Me  1  -id  i  :• 
removed.  The  Hertzian  theory  underestimates  the  contact  radii  in  each 
case  as  expected. 

The  indentation  radius  vs  impact  velocity  predict  ions  are  plotted  in 
Figures  17-19  where  they  are  compared  with  experimental  data.  In  each 
case,  the  predictions  overestimate  the  indentation  radii.  Again,  it  seems 
likely  that  the  variation  in  flow  stress  with  strain  rate  is  causing  the 
material  to  respond  more  elastically  than  it  does  during,  static  indenta- 
t  ion  . 

The  predicted  impact  load  values  are  plotted  versus  impact  velocity 
in  Figures  20-2 2.  The  loads,  calculated  using,  the  e  I  as  t  i  c- p  1  as  t  t  c  and  the 
modified  elastic-plastic  theories,  are  higher  than  the  Hertzian  theory 
values  ref  lei  ting  the  fact  that  the  contact  radii  are  overestimated  as 
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shown  in  Figures  14-16  so  that  when  the  contact  pressures  are  summed  over 
the  contact  area,  the  loads  are  overestimated.  As  in  the  case  of  ZnS,  one 
way  to  obtain  more  realistic  estimates  is  to  use  the  experimental  (impact) 
contact  radii  together  with  the  experimental  contact  characteristics  k  and 
£  to  estimate  the  load.  However,  unlike  the  case  of  ZnS ,  when  this  was 
done  the  loads  calculated  using  the  elastic-plastic  theory  did  not  show  the 
expected  reduction. 

As  mentioned  previously,  the  coefficient  of  restitution  responds 
very  sensitively  to  the  variations  in  the  energy  absorption  mechanisms. 

In  Figures  23-25,  predicted  values  of  coefficient  of  restitution  at  various 

impact  velocities  are  compared  with  experimental  data.  Again,  the  modified 

theory  yields  the  best  results  mainly  because  in  this  case  the  k/V.'  ratio 

is  one  whereas  in  the  original  theory  <l'f  1  >  1  and  the  results  are  very 

sensitive  to  this  value. 

The  radial  crack  lengths  were  measured  in  the  surfaces  of  the  stat¬ 
ically  loaded  TT-ZrO specimens,  from  the  edge  of  the  indentation  to  the 

tip  of  the  longest  crack.  The  results  are  plotted  as  suggested  by 
1/2 

P/d^  <*  c  in  Figures  26  and  27.  The  data  for  the  as  fired  TT-ZrO^, 
plotted  in  Figure  26,  show  a  linear  variation  and  the  data  points  for 
c  >  d .  fit  the  line 

i 

5  1  /2 

P/d.  =  7.93  x  10  +  1.41  r  (11) 

l 

with  a  coefficient  of  determination  of  0.695.  The  data  points  for  c  •  d . 
were  not  used  for  the  computation  because  near  the  contacts  the  crack 
lengths  are  considered  to  he  transitional  and  mav  varv  because  el  other 
charac t  er i s  t  i cs 


of  the  eon  t  ac  t  s . 
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The  line  does  not  extrapol at e  through  the  origin  hut  instead  inter¬ 
cepts  the  vertical  axis  at  an  appreciable  load.  In  chic  v; "1  1  h'  r 

TT-ZrO.j  is  somewhat  different  from  the  XnS.  This  dilference  has  not  yet 
been  explained. 

The  data  for  the  as-fired,  under-aged  and  peak-aged  specimens  are  com¬ 
pared  in  Figure  27.  Although  the  data  are  somewhat  sparse  and  scattered, 
the  results  fail  to  show  the  expected  decrease  in  crack  length  wiLh 

increasing  Kr  that  occurs  as  a  result  of  aging.  The  K,  values  ranging 
J.  c  1  c 

1/2 

from  4.5  to  6.5  Ml’am  were  expected  to  yield  a  two  to  one  variation  in 
crack  lengths. 

The  failure  to  benefit  from  increased  KT  with  aging  in  TT-ZrO.  was 

1  c  l 

reported  previously  .  In  that  case  it  was  found  that  the  number  of 
radial  cracks  decreased  with  increasing  aging.  New  data  supporting  this 
observation  are  presented  in  Figure  28.  it  is  well  known  that  the  stress 
intensity  factor  at  a  group  of  radiating  cracks  decreases  with  increasing 
number  of  cracks  (decreasing  spacing).  Therefore,  it  is  reasonable  to 
expect  that  the  increased  number  of  cracks  in  the  materials  with  lower 
might  serve  to  limit  the  crack  propagation  compared  with  that  otherwise 
expected,  yielding  results  like  those  shown  in  Figure  27. 

Another  possible  explanation  is  that  residual  stresses  increase  radial 
crack  propagation  and  that  the  residua]  stresses  increase  with  increased 
aging,  thus  offsetting  the  effect  of  the  increased  .  This  argument  is 
contradicted  by  the  fact  that  there  is  no  appreciable  variation  in  plastic 
flow  properties  with  aging  (Figure  13). 

t’sing  the  statically  determined  crack  length  relation  ( lupin  t  ion  11) 
together  with  Srpial  ion  (10),  the  theoretical  crack  loncth  vs 
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curves  were  plotted  and  compared  with  the  experimental  data  in  Figure  29. 


The  original  elastic-plastic  theory  drastically  overestimates  the  crack 
lengths.  However,  the  modified  theory  yields  much  more  reasonable  results, 
slightly  overestimating  the  crack  lengths  at  low  velocities  and  under¬ 
estimating  the  crack  length  at  the  highest  velocity. 


Radial  Crack  Length  -/xm 
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VI.  Summary  and  Discussion 


The  analysis  developed  in  this  program  represents  an  attempt  to 
obtain  a  general,  self  consistent  means  to  calculate  a  wide  range  of 
impact  characteristics  using  easily  measured  contact  characteristics 
obtained  under  static  loading  conditions.  Empirical  equations  are  used 
to  relate  the  contact  and  indentation  pressures  to  the  normalized  contact 
and  indentation  radii.  These  empirical  equations  are  manipulated  using 
Newtonian  menhuni  ,  Hertzian  '•!  ast  i  ci  ty  and  Traci  ur“  mechanics  method:-, 
to  obtain  equations  for  the  contact  radius,  indentation  radius,  load, 
penetration  time,  crack  size,  remaining  strength,  contact  time,  and 
coefficient  of  restitution. 

The  comparisons  of  the  theoretical  calculations  and  experimental 

data  presented  in  this  report  show  that  the  results  of  the  analysis  can 

be  used  to  make  reasonable  predictions  of  the  response  of  a  variety  of 

ceramic  materials  to  low  velocity  impact  by  spheres.  On  the  whole,  the 

results  were  improved  by  substitution  of  K '  for  k  in  the  various  equations 

to  obtain  the  so-called  modified  elastic-plastic  theory.  further  evidence 

1/2 


of  the  applicability  of  the  relation  I’/z  "  c 

o 


to  the  localized  impact 


damage  case  was  presented.  The  exact  form  of  the  relation  can  have  an 

important  effect  on  the  results  when  crack  lengths  are  predicted  hv  ext rn- 

3/2 

polation  to  higher  velocities  where  the  commonly  used  P  ,T  c  relation 
tends  to  underestimate  the  crack  length. 

The  f ractograph i c  investigation  reported  recent  and  the  present 

report  have  revealed  a  number  of  cases  in  which  the  present  mathematical 


model  falls  to  provide  a  completely  satisfactory  representation  ot  the 
impact  phenomena: 

1.  The  extent  of  plastic  deformation  seems  Lo  depend  on  the  strain 
rate  dependence  of  the  flow  stress.  Therefore,  the  contact 
characteristics  determined  by  static  indentation  provide  only 

an  approximate  representation  of  the  plastic  deformation  occurring 
during  impact.  The  differences  in  extent  of  plastic  deformation 
have  important  effects  in  the  partitioning  of  energy  among  the 
various  energy  loss  mechanisms,  and  the  absolute  and  relative 
sizes  of  the  various  types  of  cracks. 

2.  The  analysis  neglects  the  stress  wave  energv  and  the  fracture 

(2  U) 

energy  contributions  to  the  energy  balance.  Although  Hutchings 
has  determined  that  the  stress  wave  energy  makes  only  a  minor 
contribution  to  the  energy  losses  during  impacts  in  which  plastic 
deformation  occurs,  it  is  a  fact  that  the  above-mentioned  strain 
rate  dependence  of  the  flow  stress  favors  responses  that  are  more 
elastic  than  otherwise  expected  and,  in  these  eases,  because  the 
overall  losses  are  lower,  the  stress  wave  energy  is  proport i onat el v 
more  important.  Kirchner  and  Cruver^^  showed  previously  that  the 
fracture  surface  energy  makes  a  negligible  contribution  to  the 
energy  Losses  for  low  velocity  impacts  of  glass  spheres  on  glass 
plates.  However,  at  higher  velocities  crushing  occurs  at  the 
impact  site  and  there  is  a  substantial  reduction  in  the  cool  I  irient 
of  restitution  indicating  a  substantial  increase  in  energv  losses. 
Under  these  conditions,  little  is  known  about  the  proportion! n-  of 
the  energy  among  the  indentation,  stress  wave,  friction,  and 
fracture  energy  loss  mechanisms. 

(2b) 

3.  Marshall  and  Lawn  have  shown  that  residual  stresses  induced 

during  static  indentations  have  a  substantial  influence  on  the 
extent  of  crack  propagation  and  remaining  strength  of  glass. 

T.  irchner  and  Isaacson have  shown  that  residua!  stresses 


induced  during  single  point  machining  of  silicon  nitride  ceramics 
have  similar  effects.  It  is  reasonable  to  expect  that  residual 
stresses  are  induced  during  localized  impact  hut  the  effects  of 
these  stresses  are  not  accounted  for  in  the  analysis. 

4.  The  differences  in  induced  by  aging  treatments  in  TT-ZrO?  do 

not  result  in  the  approximately  two  to  one  differences  in  crack 

size  expected  based  on  the  present  analysis.  Failure  to  observe 

these  differences  may  occur  because  the  increased  numbers  of  cracks 

induced  in  the  lower  KT  materials  mav  act  to  reduce  the  K,  at  the 

Tc  1 

tips  of  these  cracks,  thus  reducing  the  crack  propagation  below 

that  otherwise  expected.  This  observation  confirms  an  earlier 

similar  observation  and  may  he  important  because  increasing  the 

Kr  has  been  recommended  bv  several  investigators  as  a  means  to 
1  c 

improve  the  resistance  to  localized  impact  damage. 

Because  the  important  role  of  the  strain  rate  dependence  of  the  flow 
stress  is  not  accounted  for  in  the  contact  characteristics,  x,  f,  and 

C',  determined  by  static  loading,  it  is  highly  important  to  develop  a 
method  of  determining  the  contact  characterist ics  by  impact  tests.  The 
present  investigation  indicates  that  the  theoretical  calculations  of  the 
coefficient  of  restitution  are  very  sensitive  to  the  •  to  •'  rat  in.  There¬ 
fore,  it  seems  reasonable  to  attempt  to  determine  the  contact  characteristics 

by  fitting  Lite  theoretical  curves  t *  1 1  ■  •  i  e  i  •  a '  r-  i  *  u‘  ;  <•»,  Is:-.. 

This  procedure  would  also  account  tor  the  stress  wave  and  fracture  energy 
losses.  Development  ol  this  suggested  method  of  determining  the  contact 
characteristics  is  recommended.  The  results  can  he  evaluated  based  on  the 
self  consistency  ol  the  calculated  values  of  the  other  impart  eharaeteris- 
t  ics. 
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